Detecting the quantum discord of an unknown state by a single observable 
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We propose a single observable to witness the nonzero quantum discord of an unknown quantum 
state provided that we have four copies of the state. The expectation value of this observable 
provides a necessary and sufficient condition for nonzero quantum discord in 2 x N systems and a 
necessary condition in higher finite-dimensional bipartite systems. Furthermore, a nontrivial lower 
bound of quantum discord can be obtained from this expectation value. The proposed observable can 
be experimentally measured in exactly the same easy manner as that of the entanglement witness, 
moreover a quantum circuit is designed to determine the expectation value of our observable with 
four simultaneous local qubit-measurements. 
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I. INTRODUCTION 

Entanglement is an essential resource in almost all 
quantum computing and informational processing tasks. 
In certain quantum computing tasks, however, there are 
still quantum advantages with the absence of entangle- 
ment. One typical example is the deterministic quan- 
tum computation with one qubit (DQC1) in which 
the quantum discord, introduced by Ollivier and Zurek 
and independently by Henderson and Vedral [|[, is 
proposed to be responsible for the quantum speedup [1] . 
Ever since, the quantum discord, a nonnegative number 
quantifying all the quantum correlations in a state be- 
yond the entanglement, has attracted much interest in 
quantum information theory (5WI3I]. Subsequently the 
operational meaning of the quantum discord has estab- 
lished firmly its status as an essential resource 

Quantum discord provides a measure for the quantum 
correlation beyond entanglement, which is defined as 

D A (g) = nun'Y]p k S(Q B \k) + S(g A ) - S(g), (1) 
{Ek} k 

where S(g) = — Tr(glog 2 g) is the von Neumann entropy, 
QA {qb) is the reduced density matrix of subsystem A 
(B), and the minimum is taken over all possible positive 
operator- valued measures {E k } on subsystem A [3( with 
p k = Tr{E k ® Ig) and g B]k = Tr A (E k ® lg)/p k . The 
minimum can also be taken over all the von Neumann 
measurements @, and these two definitions coincide in 
the case of zero quantum discord. In general, the quan- 
tum discord is asymmetric, i.e., D A (g) is not equal to 
Db{q) which is obtained by measuring subsystem B. For 
simplicity we only consider D A {g) in the following. 

Though almost all quantum states have nonzero quan- 
tum discord 15), states with vanishing quantum discord 
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find numerous applications in fundamental concepts in 
addition to its initial motivation in pointer states 0. 
For examples, the vanishing quantum discord is related 
to the complete positivity of a map [16( and the local 
broadcasting of the quantum correlations [l7], HI] ■ For a 
given bipartite state in a finite-dimensional system, suf- 
ficient conditions as well as necessary and sufficient con- 
ditions for nonzero quantum discord have been proposed 
[l9l - l22| in the form of local commutativity, strong posi- 
tive partial transpose, etc. Specially Rahimi and SaiToh 
[23j have introduced an example of nonlinear witness for 
the nonzero quantum discord. All these criteria, however, 
require some knowledge of the quantum states and have 
nothing to say about the nonzero values of the quantum 
discord. The problem of witnessing the nonzero quantum 
discord of an unknown state remains open. 

In this paper we shall propose a single observable to 
detect the nonzero quantum discord of a completely un- 
known state provided that we have four copies of the 
state. The expectation value of this observable provides 
not only a necessary and sufficient condition for nonzero 
quantum discord in 2 x N systems and a necessary con- 
dition in higher finite-dimensional bipartite systems, but 
also quantitative lower bounds for the quantum discord 
and the geometric measure of quantum discord. Our ob- 
servable can be experimentally measured in exactly the 
same easy manner as that of the entanglement witness, 
and moreover a quantum circuit is designed so that its 
expectation value can be easily determined via four si- 
multaneous local qubit-measurements. 



II. SINGLE-OBSERVABLE DETECTION OF 
QUANTUM DISCORD 

It turns out that, similarly to entanglement witnesses, 
there is a Hermitian observable W that detects perfectly 
the nonzero quantum discord without state tomography 
and provides a lower bound for the quantum discord. To 
construct this observable we note firstly an important 
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lesson learned from 23[ that a linear witness can never 
detect nonzero discord of a separable state. Therefore 
nonlinear witnesses should be used to detect quantum 
discord and the expectation values of an observable on 
multi copies of a state is naturally nonlinear. Secondly, 
since the quantum discord is invariant under local uni- 
tary (LU) transformations, the observable witnessing the 
quantum discord should also be LU-invariant. 

In general, the polynomial LU invariants of degree 
k for a bipartite given quantum state qab is given by 
Tr(Ug® k ), where U is a tensor product of some permu- 
tation operator acting on copies of subsystem A and some 
permutation operator of subsystem B [24| . For later use 
we introduce the following four LU invariants Tr([/jp® 4 ) 
(z = 1,2,3, 4) of order 4 where 



Us = v£v£v£V&, 



u 2 = vdvgvil, 

U4 = VflT^f F 3 f , 



(2) 



in which we have denoted by V i 



A.B 



H kl \kl)(lk\f- B the 



swapping operator acting on the i-th. and the j-th copies 
of subsystem A (B). Our discord witness reads 



where fio denotes the set of zero-discord states and \\X — 
Y\\ 2 = Tr[(X — Y) 2 ] is the square norm in the Hilbert- 
Schmidt space. 

Remark 3. Suppose that {Ai}^ and {Bj}j^ are 
complete sets of local orthogonal observables (LOOs) 
for subsystem A and B (25[, respectively, which satisfy 
Ti(AiAj) _ 

B Q = l B /^/a ^ B. Denote that r l} = (A, 
R := (rij) for i = 1, • • • , d\ — 1 and j = 0, • 



Tr(BiBj) = 5ij with A = l A / y /d A ~ and 



a {d 2 A — 1) x {d 2 B ) dimensional matrix. Thus, 



>4 - 



and 
1 is 



E E 

3=0 



Bi 



dA 



QB- 



(8) 



It can be obtained that Tr(RR T ) = Eiii Ej=o" i A i ® 
Bj) e (Ai ® Bj) e = Tr^Fj q® 2 ), where we have used 

V B = E^lo" 1 B 3 ® Bj and defined V A := eSi" 1 A ® 
Aj = U A -l A ®l A /(i A . Similarly, we have Tr[(i?i? T ) 2 ] = 

Eft! EitoVi ® Bj)c<^* ® ^)e(^* ® ® 
Bi)„ = TriVdV^VgV&Q®*). Therefore, 



W = Ui - U 3 



dA 



(U 2 - Ui), 



(3) Tr[(RR T ) 2 } - [Tt(RR t )] 2 = Tr{W g m ) = Tr{Wg m ),(9) 



where dA is the dimension of Ha and obviously 
Tr(IUp® 4 ) is also an LU-invariant. We are now ready 
to present our main result with the proof shown in the 
appendix: 

Theorem. For an unknown bipartite state g of a 
dA x ds system with 4 copies, if Tr(IU0® 4 ) = then 
D A (g) = 0, conversely if D A (g) = then Tr(Wg® 4 ) = 
when dA = 2. Moreover, if one restricts to von Neumann 
measurements in Eq. (JlJ, 

D A {g) > max{0, S{g A ) - S(g) - \og 2 [d A Q(g)}} (4) 

can be obtained, where 



Q(q) 



d A 



+ (d A ~ l)\/Tr(IU^ 4 ) + (Tre 2 ) 2 (5) 



with q := q — \aI dA ® gB- 

Remark 1. The expectation value of W provides not 
only a necessary and sufficient condition for nonzero 
quantum discord in 2 X N systems, but also a necessary 
condition in M x N systems with M > 2. 

Remark 2. Similarly to Eq. (j4)), one can also present 
a nontrivial lower bound for the geometric measure of 
quantum discord using the expectation value of W, i.e., 



Df{g) > max{0,Tr e 2 - Q(g)}. 



(6) 



D A '(g) is the geometric measure of quantum discord de- 
fined as m 



D 



A ] (e) = mm ||e-xll 



(7) 



where W = - V&V&V&VJI and the sec- 

ond equality holds since the invariance of g® under arbi- 
trary permutation of different copies. From Eq. ([9]), one 
can see that W is related to the rank of R and the ex- 
pectation value of W is always non-positive for arbitrary 
states. 

Let us take the output state of DQC1 model as an 
example. The model of DQC1 introduced by Knill and 
Laflamme can evaluate the normalized trace of a unitary 
matrix efficiently [TJ. Initially there is a single control 
qubit in the state (l+cter 2 )/2 with a describing the purity 
of the control qubit and a collection of n qubits in the 
completely mixed state l n /2 n . After the Hadamard gate 
applied to the control qubit and a control n-qubit unitary 
gate U n , the output state becomes 



gout 



2 n+l 



In 

aU„ 



aUt 

In 



(10) 



By making the measurements in the eigenbasis of o~ x and 
o~y on the control qubit, the quantum circuit provides 
the normalized trace of U n , r — Tr(C/ ra )/2 n , since the 
expectation values of o~ x and o~ y give estimates for the 
real and imaginary parts of r, respectively. 

There is strong evidence that DQC1 task cannot be 
simulated efficiently using classical computation (26| , and 
the quantum discord has been proposed as a figure of 
merit for characterizing the resources present in this com- 
putational model [ij . It is worth noticing that the output 
state (fTOf is a natural 2 x N state under the bipartite split 
between the control qubit and the collection of n qubits. 
Therefore, we can derive an explicit condition for char- 
acterizing the correlations in the output state using the 
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FIG. 1: Quantum circuit to measure the quantum discord 
witness W. Four auxiliary qubits, initially prepared in state 
|+), are control qubits for controlled-[/i (i = 1,2,3,4) oper- 
ations. Each Ui is composed of several controlled-swapping 
gates which are represented by connected crosses. After all 
the controlled-t/i operations, we perform cr^-measurements on 
the auxiliary qubits and via Eq. (|13|) in the text value of 
Tt{Wg m ) can be obtained. 



FIG. 2: Illustration of direct measurement of the observable 
W using four copies of g. There are two measurement settings 
PA4P-23P-12P-34 and P* 12 P- 3i P-i2P-34 corresponding to 
(a) U\ and U2, and (b) Us and U4. 



theorem. After some algebra, we have 



22n+3 



1, • • • ,4) operations with each qubit as source, we pre- 
form c^-measurements on all four auxiliary qubits so 
that four expectation values are obtained simultaneously. 
As it turns out (cr*) = Tr(J7i£>® 4 ), (o*) = Tr(U 2 g m ) 7 
(4) = i(Tr(C/ 3 ^ 4 ) + and (0*) = ±(Tr([/ 4 ^ 4 ) + 

(c 2 ,)), from which one can easily calculate the value of 
Tr(Wg m ) according to 



According to our theorem, DA{g ut) — if and only if 
Tr(We®? t ) = 0, i.e., a = or U% = e^l„ for some real 
(f>, which is equivalent to the result shown in Ref. [20]. 

On the other hand, if TrU n = and a = 1, using Eq. 
(U) we have a lower bound 



D A (g out ) > 1-loga I 1 + ^2 2 "+ 2 Tr(^r t ) + 1, (12) 

for the quantum discord in DQC1. It is obvious that the 
negative expectation value Tr(W gf^t) provides a non- 
trivial lower bound for DA(g ut)- It has been shown that 
for almost every unitary matrix U n (random unitary) the 
discord in the output state (p~0|) is nonvanishing 4] . How- 
ever, in the nontrivial case C/ 2 = e z ^l„ the quantum dis- 
cord DA{g ut) = 0, which provides a counter example 
for the nonvanishing quantum discord of every nontrivial 
unitary matrix U n . Thus, the quantum discord might 
not be the reason for the quantum speedup in DQC1. 



III. MEASURING THE DISCORD WITNESS W 

In order to measure the expectation value of the 
observable W we have designed a quantum circuit as 
shown in Fig. [TJ We have introduced four auxil- 
iary qubits, which are all initially prepared in the state 
|+> = (|0) + |l))/\/2. After the controlled- U, (i = 



Besides the quantum circuit shown above, one can also 
directly measure the single observable without the quan- 
tum circuit. Let P_ = (1 — V)/2 be the antisymmetric 
projector on some two copies and we have V = 1 — 2P_ 
so that we can expand W with P—ij as what follows 



W 



' i i — 2P_ 14 ) 



(8(1^2 -12) 

-(lf 2 -2Pf 12 ) 



d.,\ 



23 



2P^ 



1P B 



(lf 4 -2P* 4 ) 



d,A 



34 



2P: 



31 



v lf 2 2P B 12 ) 



(lf 4 -2P* 4 ). 



(14) 



Since the 

pA pA pB pB 
r -14 r -23 r -12 r -34 

observables like P^P^^P®^ 



measurement of the observable 
gives also information about 
etc., we need only 
two measurement settings P\^P^ 2 ?,P~i2P~d A an d 



pA pA pB pB 
r - 1 2 r - 34 r - 1 2 r - 3 4 



23- 1 -\2 ± -34 

to measure the observable W. We 
illustrate these two measurement settings in Fig. [21a) 
and Fig. [lib), respectively. 

For qubits, P_ takes the particularly simple form 
P_ = \ip-){i)-\ with \%l>~) = (|01) - |10))/V2, and it 
has already been measured in several experiments using 
two copies of g. For example, there are at least two meth- 
ods to measure P_ in photonic system. The first method 
to project two photons onto the singlet state is to use 
a Hong-Ou-Mandel interferometer 27[, which has been 
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used in (28|. The second method, employed in (29[, is to 
distinguish the Bell states with a controlled-NOT gate, 
which can transform the Bell states to separable states. 

Each of the two approaches proposed above to obtain 
the value of Tr(W g® 4 ), namely by the quantum circuit 
and the direct measurement, has its own advantages. The 
quantum circuit requires only one local measurement set- 
ting, but it needs the help of four auxiliary qubits and 
four controlled-C/ gates. The direct measurement of W 
does not require any auxiliary qubit, however, it needs 
two measurement settings. Compared with the quantum 
circuit, the direct measurement of W might be a little 
easier to realize in experiments, since the measurement 
of P- for qubits has been realized in experiments. 

IV. DISCUSSION AND CONCLUSION 

In this work, we present a single observable W to 
witness the quantum discord of an unknown state pro- 
vided there are four copies of the state. The expecta- 
tion value of our observable provides a necessary condi- 
tion for nonzero quantum discord which turns out to be 
also sufficient in 2 x JV systems. Also a nontrivial lower 
bound, which is illustrated in the example of DQC1, is 
placed by the discord witness. Moreover, a quantum cir- 
cuit is designed to determine the expectation value of the 
discord witness W using four simultaneous local qubit- 
measurements. Interestingly, four copies of states have 



also been turned out to be sufficient for entanglement 
detection of two-qubit states [30j. 
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APPENDIX 

We will present the calculation details to support our 
results and statements in the main text. 

Complete proof of the Theorem. Suppose that 
{j4j} i= f 1 and {Bj}jjL 1 are complete sets of LOOs for 
subsystem A and B (25[, respectively, which satisfy 
Tr (Ai A j) = Tr(BiBj) = with A = l A /y/d~A~ and 
Bq = 1 B l^fdji- As an example, one can choose A§ — 

1 A /Vdl, Bo = 1 B /Vd^, and {B^fg' 1 as 

the normalized traceless Hermitian generators of SU(cU) 

{^i\i=i and SU(d B ) {Ajj^j -1 , respectively, i.e., TrA 4 = 
TrAj = and Tr(XiXj) = Tr(AjAj) = Therefore any 
bipartite density matrix has expansion 



1 A 1 B V 
-r- ® -r- + > 

d A d B ^ 



Xj\i 



+ 



.7 = 1 



Vj 



d A -ld B -l 

E E 

i=l .7=1 



(15) 



where the coefficients Xi, yj, tij forms the column vectors 
x, y and the {d\ — 1) X (d 2 B — 1) dimensional matrix T, 
respectively. 

Denote by R a (d\ — 1) x (d B ) dimensional matrix 
whose matrix elements are r,- 



I,-- - ,d% - 1 and j 
given in Eq. (|15p . we have 



0,--. ,d% 



(Ai (g> Bj) e with i 
1. For the state g 



R= (x (d 2_ 1)xl | r (d 2_ 1)x(d |_ 1) 



In general we have 



(16) 



where roj '■= (Aq ® Bj) B and we have used 



4-i 

^2 r ojA <g) Bj 

3=0 

E r j^==®Bi 

j=o 
4-i 



4-i4-i 
Q = E E n 3 Al ® B 3 

i=0 j=0 

4-i4-i lA 

= E r ijA i ®B j + — ®QB, (17) 

i=i j=o A since r oj 



— i 

dA 



idA 

QB 



(18) 



{1 A /Vdl ® Bj) e = Tr(g B B 3 )/y/dA~. A 
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straightforward calculation yields 
Tr(RR T ) 

= e r io r ij 

4-i 4-i 

1=1 3=0 



(19) 



where A™ and -B™ operate on the mth subsystem A and 
the nth subsystem B, respectively. We have used the 
swap operator 



Vii= £ B > B . 

3=0 



(20) 



which operates on the first and second copies of subsys- 
tem B, and defined 



Therefore, 



TtKRR 1 Y] - [Tt(RR )Y = Tv{W'g m ) = Tr(Wg m p6) 



where we have defined 



and the second equality in Eq. f[26]) holds since 

Tr(VdV* V£ g^) = Tr(V^f F 3 f f?® 4 ), (28) 
^(V&V&V&q**) = Tr(V 3 ivBv£ g ®% (29) 

which comes from the invariance of g® 4 under arbitrary 
permutation of different copies. 

(i). If Tv(Wg m ) = then we have Tr[(RR T ) 2 } = 
[Tr(RR T )} 2 , which means that the rank of the matrix R 
is or 1. 

(ia) If rank(i?) = 0, then R = 0. From Eq. |jnj>, we 
can obtain 



(li-i 



V 



V 



d A 



(21) 



operating on the first and second copies of subsystem A. 
Notice that the swap operator is defined as 



7:=£|fc)<*|®|l><fc|. 



(22) 



A7 



It is worth noticing that every Hcrmitian operator can 
be written in the operator basis, for example, the swap 
operator V\f on the first and second copies of subsystem 
B can be written as 



4-i 
= £ ClJ B} ® B 



(23) 



where 



cy = TrOB^BfVii) 
= ^{ll^lfcXfcl^ll) 



(24) 



Q d A 



Qb, 



(30) 



which obviously has zero quantum discord D A (g) = 0. 
(ib) If rank(i?) = 1, then there exist real numbers a,; 



and j3j such that r. 



with i = 1, 



j = 0, • • • ,dr B -l. It follows from Eq. JI7| that 



^ — 1 and 



4- 1 
£> = E 



4-i 

E h B i 

3=0 



it 

d A 



(31) 



Using the eigenvectors {11^} of J2i=i 1 a iAi one can get 

£ fc (n fc ® i)e(n fc ® l) = & i.e., - o. 

Therefore, if Tr(IU£>® 4 ) = then we have D A (g) = 0. 

(ii). Conversely if d A = 2 and D A (g) = then there 
exists an ortho normal basis {IV'i), IV^)} of subsystem A 
and two density matrices g k of subsystem B such that 

g = p\ih)(ih\ ® ef + (i -p)I^X^I ® f?f , (32) 



Hence, one can directly obtain Eq. (|20|) using Eqs. ([23 
and (|24p . On the other hand, we have 

Tr[(i?i? T ) 2 ] 

4-i4-i 

= E J2( A i® B i)e( A l® B ])e( A l® B l)Mt® B , 

i,k=l j,l=0 

= T^iviviV^). (25) 



with < p < 1. Using Eq. (f3"2"j). it is straightforward to 
check that Tr(IUp® 4 ) = 0. Therefore, for a 2 x N state 
ft the discord D A {g) = if and only if Tr(Wg m ) = 0. 

Moreover, if we perform a von Neumann measurement 
{IIfc} fc f. 1 on the subsystem A then we will obtain the 
outcome k with probability p k = Trfllfc® Ig) leaving the 
subsystem £? in the state ftg|j, = g' B , k /pk where g' B , k = 



6 



Tr a (11^ (g> lg). The conditional entropy becomes 



d.A 

k=l 
d A 

k=l 

d.A 

= -J2 T ^ gl B\k l °S2(0'B\k/Pk)} 
k=l 

d.A d A 

= J2Pklog 2 Pk-J2Tr{g' Blk log 2 g' Blk ) 

k=l k=l 
d A 

= E Pk lo &2 Pk + s (Efcii n fc ® e' B 



k=l 



> -log 2 ^-log 2 {EfciiTr[(e' B|/£ ) 2 ]}, (33) 

where we have used the inequalities Efcii Pk ^°§>2Pk > 
— log 2 dyi and S(g) > — log 2 Trg 2 . By expanding Ilfc 
using a complete set of LOOs with coefficients e\ = 

Tr{H k Ai), i.e., U k = £<f _1 $A u and denoting by M 
a (d\ — I) x (d 2 A — I) matrix with elements My = 
Tr^ 1 ® AfVj p® 2 ), we obtain 



dA 

Y, T ^B\k) 2 ] 

k=l 

dA 

E T ^'Jr* 



k=l 

fc=l i ,j=o 

dA Trn 2 



fc— 1 ij — l 
d A d 2 A -l 



dA 



< EE ^V^mm^ 



k=l i=l 



(dA - l)v/Tr(W^® 4 ) + (Trg 2 ) 2 + ^ 

«A 

Q(e), 



(34) 



where we have used the Cauchy-Schwarz inequality, 



" (ef) 2 = - l)/d A , and Efc=r e? = for 
i = 1, • • • ,d\ — 1. Therefore, 

£U(e) 

= SW) - S{g) + min J2kPk s (QB\k) 

{ E ki 

> S(g A ) - S(g) - log 2 d A - log 2 { £^ Tr[(g' B|/£ ) 2 ]} 

> S-( eA ) - S(g) - log 2 d A - log 2 [Q(f?)] 



which is the lower bound of the quantum discord Eq. Q 
in the main text. 

Remark. Similarly to Eq. (|35|) . one can also present 
a nontrivial lower bound for the geometric measure of 
quantum discord D A (g) using the expectation value of 

W, where D A '(g) is the geom etric measure of quan tum 
discord defined in Ref. [2(|. Actually, Ref. [31j has 
shown that 



^ 2) (p) = min{Tre 2 -E fe Tr(e' B|/£ ) 2 }, 

Using Eq. (|3"4"1) . one can obtain 

Df{g) = min{Tr, 2 -E fe Tr(^ |fe ) 2 } 
= Trg 2 -max{E fc Tr( ' B|fc ) 2 } 



(36) 



> ^g 2 ~Q(g), 



(37) 



which is Eq. in the main text . 

DQC1 Model. The output state of DQC1 model can 
be written as 



Qout 



1 ( In aUl 
2»+! ^ aU n 1„ 

= 2^r(li ® In + «|1><0| ® tf„ + a|0)<l| ® Ul) 

= ^lE « (i+J> ° d2 K>01 ® CM*. (38) 

Notice that for DQC1 model d^ = 2, thus 

W = Z/j. + Z7 4 - U 2 - U 3 - (39) 
Based on Eq. (|3"8")l . one can directly calculate as follows, 
f 



Tr(U 2Q f: t ) = 
Tr(C/ 3 ^ M 4 t ) = 



i2n+2 



24n+4 V 
+2a 4 |Tr[/ 2 | 2 

24n+4 V 

— (2 2 

n+4 V Z 



-,2n+2 



24n+ 

2 : 

f 



8a 2 |Tr[/„l 2 



. 2 2 "+ V) 



8a 2 |Tr£/„| 2 ) 



2 2n+3 a 2 



2 2 "+V) 



! )- 



Tr(C/ 4 ^ 4 ) = ^ Ti (2 2 " +3 + 2 2 "+ 3 ^ > > 
Therefore, 

Tr(Wgf* t ) = Tr[(C/! + £7 4 - ET a - C/ 3 )^ 4 )] 



(40) 
(41) 

(42) 
(43) 



a 



22n+3 



Tr£/ 2 



(44) 



5(eA)-5( e )-io g2 [dAQ(e)], 



(35) 



which is Eq. (fTTj) in the main text. According to our 
theorem, D A {g out ) = if and only if Tr(W^f u 4 ) = 0, i.e., 
a = or U 2 = e^ln for some real 4>, which is equivalent 
to the result shown in Ref. 1201. It has been shown that 
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for almost every unitary matrix U n (random unitary) the 
discord in the output state (1551) is nonvanishing [1] . How- 
ever, in the nontrivial case U% = e^l n the quantum dis- 
cord Dji(g ou t) — 0, which provides a counter example 
for the nonvanishing quantum discord of every nontrivial 
unitary matrix U n . Thus, the quantum discord might 
not be the reason for the quantum speedup in DQC1. 

On the other hand, if TrJ7 n = and a = 1, then 
for the output state Eq. (|38|) one can obtain that 
qa = li/2, qb — ln/2", and the eigenvalues of g out are 
{0, • • • , 0, 1/2", • • • ,1/2"} in which there are 2" copies of 
and 2" copies of 1/2". Therefore, 

S{qa) = 1, (45) 
S(Qout) = n, (46) 

= ± (47) 
TreL = 2^T- (48) 

Using Eq. ()35|) we have a lower bound for the quantum 
discord in DQC1 

D A (g out ) > l-log 2 (l + ftwTrQVQt*) + lY (49) 

which is Eq. fj 12[) in the main text. It is obvious that the 
negative expectation value Ti(W gf ut ) provides a nontriv- 
ial lower bound for DA(g ut)- 

Measuring the discord witness W . In order to 
measure the expectation value of the observable W we 
have designed a quantum circuit as shown in Fig. Q] 
in the main text. We have introduced four auxiliary 
qubits, which are all initially prepared in the state |+) = 
(|0) + |l))/\/2. After the controlled-^ (i = ,4) 
operations with each qubit as source, we preform ad- 
measurements on all four auxiliary qubits so that four 



expectation values are obtained simultaneously. As a re- 
sult, the total state including the four auxiliary qubits 
is 

Qtotal = le E (\i}(i'\®\j)(f\®\k)(k'\®\W\ 

mlU 3 k mU[gf u \ufufufuf). (50) 
Thus, we have 

(al) = i[Tr(^ 4 C/ 1 t ) + Tr(C/ 1 ^ 4 )] 

= Tr(l^ 4 ), (51) 

(o*) = i[Tr(^ 4 [/t) + Tr(C/ 2 ^ 4 )] 

= Tr(U 2 g m ), (52) 

(4) = i[Tr(^ 4 C/ 3 t ) + Tr(C/ 3 ^ 4 ) + 2(ai)] 

= llTr(U 3 g^) + (al)}, (53) 

(a 4 ) = i[Tr(^ 4 C/|) + Tr(C/ 4 ^ 4 ) + 2(^)] 

= i[Tr(l^ 4 ) + (^)], (54) 

from which one can easily calculate the value of 
Tr(Wg m ) according to 

Tr(Wg^) = 2(Vi) - (at)) -^-({a 2 x ) - ( CT 4 >), (55) 

which is Eq. (ll3l) in the main text. In order to obtain 
Eq. (IT4|) in the main text, one can use the definition of 
W Eqs. ((2]) and (|3]) in the main text and the equation 
V = I - 2P-. 
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